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An  optimal  maneuver  strategy  considering  terminal  guidance  accuracy  for  hypersonic 
vehicle  in  dive  phase  is  investigated  in  this  paper.  First,  it  derives  the  complete  three- 
dimensional  nonlinear  coupled  motion  equation  without  any  approximations  based  on 
diving  relative  motion  relationship  directly,  and  converts  it  into  linear  decoupled  state 
space  equation  with  the  same  relative  degree  by  feedback  linearization.  Second,  the  diving 
guidance  law  is  designed  based  on  the  decoupled  equation  to  meet  the  terminal  impact 
point  and  falling  angle  constraints.  In  order  to  further  improve  the  interception  capability, 
it  constructs  maneuver  control  model  through  adding  maneuver  control  item  to  the 
guidance  law.  Then,  an  integrated  performance  index  consisting  of  maximum  line-of-sight 
angle  rate  and  minimum  energy  consumption  is  designed,  and  optimal  control  is 
employed  to  obtain  optimal  maneuver  strategy  when  the  encounter  time  is  determined 
and  undetermined,  respectively.  Furthermore,  the  performance  index  and  suboptimal 
strategy  are  reconstructed  to  deal  with  the  control  capability  constraint  and  the  serous 
influence  on  terminal  guidance  accuracy  caused  by  maneuvering  flight.  Finally,  the 
approach  is  tested  using  the  Common  Aero  Vehicle-H  model.  Simulation  results  demon¬ 
strate  that  the  proposed  strategy  can  achieve  high  precision  guidance  and  effective 
maneuver  at  the  same  time,  and  the  indices  are  also  optimized. 

©  2014  Published  by  Elsevier  Ltd.  on  behalf  of  IAA. 


1.  Introduction 

Inspired  by  the  excellent  performance  of  hypersonic 
vehicle,  such  as  flexible  ballistic  form  and  high  lateral 
maneuver  capability,  much  attention  has  been  played  on 
the  corresponding  research  in  recent  years.  As  a  popular 
research  topic  in  aerospace  field,  the  typical  representation 
includes  the  American  common  aero  vehicle  (CAV)  and 
Hypersonic  Technology  Vehicle  2  (HTV-2)  which  is  con¬ 
sidered  as  the  expansion  of  CAV  [1,2],  Dive  phase  is  the 
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last  one  in  the  whole  flight  phases,  and  the  guidance  is 
concerned  with  steering  the  vehicle  to  reach  the  desig¬ 
nated  termination  target  with  prescribed  accuracy.  Due  to 
the  highly  coupled,  rapid  varying,  nonlinear  motion  model 
which  contains  large  parametric  uncertainties,  the  gui¬ 
dance  which  can  satisfy  multiple  terminal  constraints  in 
presence  of  necessary  path  constraints  and  the  maneuver 
penetration  remain  the  key  technical  challenges  [3],  Con¬ 
sequently,  the  diving  maneuver  strategy  design  consider¬ 
ing  the  terminal  guidance  accuracy  is  more  complicated 
and  significant. 

Most  current  researches  focus  on  the  guidance  law 
derivation  and  maneuver  model  design  separately.  As  for 
the  high  precision  diving  guidance,  traditional  algorithms 
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simplify  the  three-dimensional  (3D)  diving  motion  into 
the  two-dimensional  planar  motion,  and  derive  the  corre¬ 
sponding  guidance  law  via  advance  control  theories,  such 
as  optimal  control  [4,5]  and  sliding  mode  control  [6], 
However,  despite  the  different  control  theories  were  fully 
employed,  nearly  all  these  guidance  laws  were  deprived 
from  zero  line-of-sight  (LOS)  angle  rate.  Therefore,  the 
trajectories  were  so  straight  forward  that  the  penetrative 
ability  will  be  greatly  declined  [7],  In  order  to  improve 
maneuver  performance,  maneuvering  flight,  such  as  weav¬ 
ing  and  barrel  roll  maneuver  can  enhance  penetrative 
ability  effectively  [8],  Great  many  factors  can  influence 
penetrative  ability,  in  which  the  maneuver  strategy  is 
the  most  important  one  [9],  Actually,  the  diving  maneuver 
penetration  process  is  the  mutual  combat  between 
pursuer  and  evader;  therefore,  both  pursuit-evasion 
and  differential  game  can  be  well  adopted  to  obtain  the 
optimal  maneuver  strategy,  in  which  the  perfect  move¬ 
ment  information  for  each  other  is  needed  [9-11].  As  a 
result,  the  strategy  will  be  much  different  if  the  movement 
information  of  pursuer  cannot  be  determined  completely. 
Using  the  limited  ballistic  parameters,  Shinar  gave  the  3D 
maneuver  control  law  [12],  Guo  et  al.  discussed  the 
maneuver  control  problem  in  ballistic  middle  phase  with 
the  assumption  that  the  Kalman  filter  was  used  to  predict 
the  movement  of  evader  and  the  evader  adopts  Singer 
maneuver  model  to  increase  prediction  error,  in  which 
the  Hoo  maneuver  performance  index  is  constructed 
with  minimum  energy  consumption  and  maximum  pre¬ 
diction  error  [13-16],  In  addition,  Ohlmeyer  analyzed  the 
influence  on  miss  distance  of  sinusoidal  maneuver,  and 
presented  an  analytical  expression  of  miss  distance  con¬ 
sidering  target  initial-conditions  uncertainties  [17],  How¬ 
ever,  it  is  known  that  maneuvering  flight  can  cause  serious 
influence  to  the  terminal  guidance  accuracy,  and  even  lose 
the  target  [18],  Most  current  research  about  maneuvering 
flight  focuses  on  maneuver  model  and  strategy,  ignoring 
the  influence  on  guidance  accuracy  caused  by  maneuver, 
which  means  that  it  is  necessary  to  construct  the  inte¬ 
grated  performance  index  and  design  the  novel  maneuver 
strategy  considering  guidance  accuracy. 

Motivated  by  the  aforementioned  problems,  a  novel 
diving  optimal  maneuver  strategy  is  considered  in  this 
paper,  whose  purpose  is  to  obtain  the  optimal  maneuver 
effect  along  with  the  high  guidance  accuracy  at  the  same 
time.  First,  the  complete  nonlinear  coupled  relative  motion 
equation,  as  the  basis  of  the  further  research  of  guidance 
law  and  maneuver  strategy,  is  constructed  without  any 
approximation.  Second,  with  the  help  of  feedback  lineariza¬ 
tion  and  feedback  control,  the  3D  guidance  law  which  can 
satisfy  multiple  terminal  constraints  can  be  also  obtained, 
and  it  is  convenient  to  achieve  maneuvering  flight  though 
adding  maneuver  control  item  into  the  original  guidance 
law.  Moreover,  considering  the  difference  when  the  encoun¬ 
ter  time  is  determined  and  undetermined,  it  will  establish 
the  optimal  performance  indices  to  evaluate  the  maneuver 
performance  and  design  the  maneuver  strategies  using 
optimal  control  for  each  case.  In  addition,  combining  the 
control  capability  limitation,  especially  the  terminal  gui¬ 
dance  accuracy,  the  performance  indices  under  different 
situations  will  be  reconstructed,  and  then  the  optimal 


strategy  and  guaranteed  cost  maneuver  strategy  will  be 
also  obtained.  Finally,  the  maneuver  parameters  can  be 
optimized  through  several  groups  of  simulations.  The  main 
contributions  of  this  paper  are  the  establishment  of  3D 
relative  motion  equation,  the  construction  of  maneuver 
control  equation,  the  design  of  the  maneuver  performance 
indices  and  optimal  maneuver  strategies,  the  reconstruction 
of  indices  and  strategies  that  the  terminal  guidance  accu¬ 
racy  and  control  capability  are  fully  considered. 

2.  Problem  formulation 


2.1.  Formulation  of  motion 


The  three-degree-of-freedom  (3DOF)  point-mass  dynamic 
under  ballistic  coordinate  system  is  adopted  [4],  For  hyperso¬ 
nic  vehicle,  owing  to  the  high  flight  velocity,  short  flight  time 
and  low  altitude  in  dive  phase,  compared  with  the  aerody¬ 
namic  force  and  gravity,  the  Carioles  inertial  force  and 
centrifugal  inertial  force  caused  by  the  earth's  self-rotation 
are  very  small,  and  the  accelerations  are  approximately 
10  2 ga  and  10  3g0,  respectively  [4],  Therefore,  the  non¬ 
rotating  spherical  earth  assumption  is  employed  in  the  motion 
equation  construction  and  guidance  law  design;  the  motion 
equations  in  ballistic  coordinate  system  are  given  in  the 
following  Eq.  (1): 


pv2SmCD 
=  “'"lifiP''” 
pv2SmCL  cos 
-  2mv 
 pv2SmCL  ' 


(1) 


where  the  position  coordinates  are  the  radial  distance  from 
the  center  of  the  Earth  to  the  vehicle  r,  the  longitude  X,  the 
latitude  tp\  and  the  velocity  coordinates  are  the  earth-relative 
velocity  magnitude  v,  the  velocity  slope  angle  6,  and  velocity 
azimuth  angle  a  measured  from  the  north  in  a  clockwise 
direction,  p  is  the  atmospheric  density;  m  is  the  vehicle  mass; 
Sm  is  the  reference  area;  g  =  pM/r2  is  the  gravity;  pM  is  the 
gravitational  constant.  CD  and  CL  are  the  drag  and  lift 
coefficients,  respectively,  and  both  these  aerodynamic  coeffi¬ 
cients  are  the  function  of  the  Mach  number,  altitude  and  angle 
of  attack  (sometimes  independent  of  altitude,  such  as  CAV-H 
[1]).  In  these  coefficients,  control  variables  include  angle  of 
attack  a  and  the  bank  angle  v. 


2.2.  Relative  motion  equations 

Traditional  decoupled  guidance  needs  to  divide  the 
diving  motion  into  diving  and  turning  sub-planes  with  the 
assumption  that  the  three-dimensional  diving  motion  can 
be  decoupled  into  longitudinal  and  lateral  directions  [4], 
However,  due  to  high  flight  velocity,  the  coupling  between 
longitudinal  and  lateral  motions  is  very  high,  especially 
when  large  amplitude  lateral  maneuver  is  adopted.  As  a 


:  al.  /Acta  Astronautica  104  (2014)  231-242 


233 


result,  as  the  basis  of  the  further  guidance  and  maneuver 
research,  it  is  necessary  to  establish  the  3D  coupled  relative 
motion  equations  to  describe  the  diving  motion  clearly.  To 
begin  with,  target  coordinate  system  O-xyz  with  the  origin 
on  the  target,  the  x-axis  pointing  north,  y-axis  pointing 
upward,  and  z-axis  pointing  east  is  established  and  the 
relative  motion  relationship  is  shown  in  Fig.  1 . 

In  Fig.  1,  M  and  T  refer  to  vehicle  and  target,  respec¬ 
tively,  MA  is  the  parallel  of  current  horizon.  It  is  known 
that  3D  diving  motion  can  be  described  in  longitudinal  and 
lateral  directions,  where  longitudinal  motion  can  be 
described  in  the  MAO  plane  and  lateral  channel  can  be 
determined  by  current  velocity  vector  v  and  the  compo¬ 
nent  in  the  MAO  plane  vd.  In  addition,  d  is  the  sight 
distance  between  the  vehicle  and  target,  AD  is  the  LOS 
elevation  angle,  yd  is  the  longitudinal  velocity  azimuth 
defined  as  the  included  angle  between  MA  and  vd,  t]D  is  the 
included  angle  between  vd  and  LOS,  r/T  is  the  included 
angle  between  v  and  vd,  yT  is  the  lateral  velocity  azimuth, 
and  An  is  the  LOS  azimuth.  Both  yT  and  An  are  measured 
from  the  north  in  a  clockwise  direction.  Then  the  following 
equations  hold: 

-v  cos  y,  cos  nD 

dlD  =  V  COS  r,T  Sin  r,D  (2) 

dirr  =  v  sin  r/T  cos  >iD 

in  which  d  is  the  rate  of  change  of  relative  distance. 
Differentiation  is  performed  on  the  second  and  third 
equations  of  Eq.  (2): 

f  d\D+<ti D  =  V  COS  i/T  sin  r/D  —  VflT  sin  rjT  sin  rfD +VilD  COS  rip  COS  I/D 
\  dApp+dApp  =  v  sin  i/r  cos  i/d+Vi/t  cos  i iT  cos  r/D-riDv  sin  i/T  sin  i/D 

(3) 

where  v  is  the  velocity  magnitude  differential.  Substituting 
Eq.  (2)  into  Eq.  (3)  will  yield  the  following  equations: 


(4) 


according  to  the  geometry  relationship: 

{’to  =  to + Yd 
ijt  =  Att-Yt 


then,  Eq.  (4)  can  be  rewritten  as 

{Ad  =  (y  -  2g)  Ad  -  v  S,B  w  an  ibqjh.  _  |>T)  _  d  yD 
An  =  (v~2l)2n--V  5111  "d  Si°  ’*p(AD  +  fr?)  +  arr 


(5) 


(6) 


Comparing  these  equations  with  the  decoupled  equa¬ 
tions  based  on  the  simplified  planar  motion  [5],  Eq.  (6) 
adds  the  coupling  factor  between  longitudinal  and  lateral 
motions:  therefore,  it  can  describe  relative  diving  motion 
more  accurately.  Here,  we  define  the  coupling  coefficient: 


Adt 


v  sin  t]j  sin  tjq 
d 


(7) 


It  is  evident  that  the  coupling  coefficient  ADT  increases  as 
tjD,  r/T  and  velocity  v  increase,  while  decreases  as  the  sight 
distance  d  increases.  As  a  result  the  coupling  coefficient  ADT 
will  be  close  to  zero  when  r/D,  ijt  and  v  are  small  and  the 
relative  distance  d  is  very  large,  which  means  the  coupling 
equations  can  be  the  simplified  into  decoupled  ones  as 
shown  in  Ref.  [4],  The  velocity  azimuth  rjD  and  t]T  can  be 
calculated: 


(r)T  =  At -a 

\  Id  =  arccos(d/v  cos  y-,  ) 

We  can  refer  to  Ref.  [4]  for  more  information  and  the 
computation  of  the  relative  motion  coefficients  AD,An , 
Ad,  An,  d,  d.  For  diving  guidance,  the  attack  direction  is  free, 
regardless  of  the  strict  constraints  for  terminal  impact 
point  and  falling  angle.  In  order  to  take  advantage  of 
control  method  to  resolve  guidance  problem,  we  begin 
with  the  introduction  of  state  variables,  control  variables 
and  output  variables  via  motion  Eq.  (6): 

ix  =  [Xi ,  x2]r  =  [Ad  +  k(AD  +  Of),  An]'1 
ii  =  [Ui,u2]t  =  [yd,7t]'  (9) 

h(x)  =  [hi(x),  h-zix)]1  =  [AD+k(A.D+0f),An]T 

where  Of  is  the  required  terminal  impact  angle,  k  >  0  is  the 
revise  coefficient  for  terminal  angle.  It  is  known  that  the 
guidance  objective  is  to  realize: 


!limxt(t)  =  0 
limx2(t)  =  0 


(10) 


denote  Tg0  =  -  d/d  as  the  flight  time  to  go,  then  Eq.  (6)  can 
be  rewritten  as 

!Xi  =  (v-t-7?+1)(Xi  —  k(Ao+0f))—  AdtX2+ADj-zj-+jT^'d 
*2=  (y+  7^  *2  —  71dt(Xi  —  k(AD  +  Of))  —  ADtYd  ~  rg/T 


dXD + dAD  =  ydAo  -  vf/j  sin  i/r  sin  yD  — yDd 
diYT  +  din  =  IdAn  -  vfjD  sin  i/T  sin  r/D  -  i]Td 


(11) 
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The  above  equations  can  be  expressed  as  the  standard 
form  of  nonlinear  control  system: 


'  x=f(x)+g(x)u 
y=h(x) 


(12) 


where 


f(x)  = 

g(x)  = 


—  k(2o  +  0f))—  AdtX2 
(v + '4„)  x2  —  A°r(x i  -  k(*o  +  Of)) 

i  Am  1 


(13) 


3.  3D  guidance  law  and  maneuver  control  equations 


is  nonsingular  at  the  equilibrium  state.  Therefore,  the  non¬ 
linear  guidance  system  satisfies  the  condition  in  Eq.  (14), 
which  means  it  can  be  decoupled  into  linear  equations  with 
the  same  relative  degree.  The  new  state  variables  are 


Z  =  T(X)  =  [h,  (x),  h2(x)f  =  [1D  +  k(HD  +  Of),  Arrf 
and  the  inputs  transform: 
u  =  p-,(*)[-Q(x)+«] 

where  v  =  [vy,  vz]T  is  the  new  control  variable,  and 
LfJiiCx)] 


Q(x)  = 


Lfh2(x) 


—k(Ao +&/))— AjjjX  2 

(v+ tfc)*2  ~^DT(Xl  —  ^(2 D  +  ®/)) 


(16) 

(17) 


(18) 


According  to  the  analysis  in  Section  1,  the  3D  high 
precision  coupled  guidance  law  design  is  still  a  key 
technical  challenge  for  hypersonic  vehicle.  The  diving 
guidance  problem  can  be  described  as  follows:  given  the 
initial  condition  of  the  vehicle,  design  guidance  law  to 
steer  the  vehicle  to  attack  the  ground  fixed  target  under 
the  flexible  multiple  path  constraints,  while  satisfying  the 
constraint  of  terminal  impact  angle  and  eliminating  the 
position  error.  Because  of  the  coupling  and  nonlinear 
characteristics  of  the  diving  motion  described  in  Eq.  (12), 
it  is  too  difficult  to  design  the  guidance  law  directly. 
Therefore,  we  can  employ  feedback  linearization  to  decou¬ 
ple  and  linearize  the  original  nonlinear  guidance  equa¬ 
tions,  and  then  the  design  of  the  guidance  law  can  be 
greatly  simplified.  Meanwhile,  feedback  control  can  be 
also  used  to  design  the  guidance  law.  Finally,  3D  diving 
guidance  law  is  obtained  through  substituting  the  feed¬ 
back  control  law  of  the  linear  equations  into  the  original 
nonlinear  system. 


Combine  Eqs.  (15)— (18),  we  can  get  the  linear  decoupled 
system: 


j  z  =  Av 

\y=cz 


(19) 


where  the  new  state  matrix  and  output  matrix  are 

m;  m  °i 


In  Eq.  (20),  both  A  and  C  are  in  the  Brunovsky  standard  form 
with  relevant  rank  [19],  and  it  is  obvious  that  the  linear 
control  system  is  controllable  completely.  With  the  linear 
equations  expressed  in  Eq.  (19),  the  diving  guidance  objec¬ 
tive  shown  in  Eq.  (10)  can  then  be  rewritten  as 


ilimx,(t)  =  limzi(t)  =  0 

t^tf  t-.t, 

limx2(t)  =  limz2(t)  =  0 


3.1.  3D  guidance  law 

Given  the  nonlinear  system  Eq.  (12)  and  the  corre¬ 
sponding  input  and  output  state  variables,  if  the  system 
satisfies 

(1 )  LgjLfhi(x)  =  0,(1  <  f  <  m,  1  <  j  <  m,  1  ^  k  <  rf  —  1) 

(2)  The  m  x  m  matrix: 


Then  the  feedback  control  law  of  state  space  Eq.  (20)  is 
v=-kz  (22) 

where  k  =  [kD,  kr]'  is  the  effective  navigation  ratio  for  which 
a  value  between  3  and  5  is  frequently  chosen  in  practice 
[13],  It  is  obvious  that  Eq.  (22)  is  sufficient  condition  to 
Eq.  (21),  which  means  guidance  objective  can  be  achieved 
with  the  guidance  law  (22). 


'LgiLp  'hi(x)  ■ 

LgmLj  1  h\(x) 

P(x)  = 

Lg  iLp-’fcfx)  • 

••  LgmLf2  - 1  h2(x) 

LgiUf-'hmfx)  ■ 

”  LgmLf"1  1  hm(x) 

being  nonsingular  at  x  =  x0,  then  the  system  is  said  to 
have  a  relative  degree  r = X["=  1  r„  and  it  is  linearizable  if 
sufficiently  smooth  function  y  =  h(x)  exists  so  that  the 
relative  degree  of  the  system  is  r  =  Z^=  ^ri  =  n  [19], 

The  relative  degree  of  system  Eq.  (12)  is  r=  2,  with 
ri  =  r2  =  1,  and  the  decoupled  matrix: 


r^,h,(x) 

(X>=  k,h2(x) 


ig2hi(x)l 
Lg2h2(x)  J 


(15) 


3.2.  Maneuver  control  equations 

Owing  to  the  complicated  outside  aerodynamic  force, 
gravity  and  the  influence  of  zero  LOS  angle  rate,  the  diving 
ballistic  is  much  different  from  the  two-body  orbit 
of  traditional  ballistic  missile  in  middle  phase,  so  it  is 
impossible  to  establish  the  maneuver  control  equations  via 
the  relative  position  and  velocity  directly.  Guidance  law 
(22)  can  steer  vehicle  to  attack  the  ground  fixed  target 
with  high  precision,  but  the  ballistic  will  be  too  straight  to 
realize  effective  perception  because  of  zero  LOS  angle  rate. 
Considering  that  the  guidance  law  (22)  takes  LOS  angle 
and  the  varying  rate  as  the  main  parameters  to  describe 
relative  motion,  the  maneuver  control  equations  con¬ 
struction  will  be  much  easier  if  the  same  parameters 
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are  adopted.  Therefore,  in  order  to  improve  maneuver 
performance,  we  can  add  maneuver  control  law  to 
Eq.  (22): 

v=-kz+um  (23) 

in  which  tim  =  [umD,  umr]T  is  the  maneuver  control  law,  so 
this  novel  control  law  can  be  then  called  integrated  control 
law.  With  Eq.  (23),  it  is  very  convenient  to  obtain  different 
maneuver  trajectories  by  designing  maneuver  control  law 
um.  The  practical  LOS  angle  rate  can  be  expressed  as 
follows  by  substituting  Eq.  (23)  into  state  space  Eq.  (19): 
z=-kz+um  (24) 


It  can  be  seen  from  Eq.  (24)  that  different  maneuver 
forms  can  be  obtained  though  modifying  maneuver  con¬ 
trol  law  um.  Substituting  Eq.  (23)  into  the  inputs  transform 
(17),  the  following  can  be  obtained: 


ii  =  p-V)[-Q(x)+o] 

>3  Jfi 

(\/T2g0)-A2m  I  -Adt 


-  (p+j^+'l^D+ADTx2-kDXi+umD 

-  (!+t^)*2  +ADTkD  -  kTx2  +  umT 


(25) 


Eq.  (25)  shows  that  the  diving  guidance  law  which  can 
realize  high  precision  guidance  and  maneuvering  flight  is 
obtained  by  substituting  the  control  law  derived  from 
feedback  control  and  additional  maneuver  item  into  the 
original  nonlinear  equations.  Based  on  the  velocity  azi¬ 
muth  rates  studied  above,  the  required  control  overload 
can  be  expressed  as 


I  ny  =  /Dv/g0 
l  nz  =  rTv/g0 


where  ny  and  nz  are  the  required  longitudinal  and  lateral 
control  overloads,  respectively,  g0  is  the  gravity  accelera¬ 
tion  on  the  sea  level.  Then  the  total  overload  is 


The  angle  of  attack  can  be  obtained  from  (28)  by 
inverse  interpolation. 

Y  =  CLqS=g0n2  (28) 

where  Y  is  the  aerodynamic  lift,  q=pv2/ 2  is  the  dynamic 
pressure  and  the  bank  angle  can  be  expressed  as 
v  =  arctan  (nz/ny)  (29) 


4.  Optimal  maneuver  strategy  when  encounter  time  is 
determined 

It  is  clear  that  the  diving  maneuver  penetration  is  the 
confrontation  with  air  defense  system.  The  defense  system 
consists  of  detective  radar  and  interceptor.  Detective  radar 
is  adopted  to  measure  and  estimate  the  movement  infor¬ 
mation  of  the  vehicle  and  transit  them  to  the  ground- 


based  defense  system,  which  can  be  used  to  steer  the 
interceptor  to  the  target.  As  discussed  in  the  above  section, 
maneuvering  flight  form  and  the  maneuver  strategy  play 
the  key  role  in  the  improvement  of  interception  capability 
effectively.  However,  we  cannot  make  certain  whether  the 
vehicle  acknowledges  the  movement  information  of  inter¬ 
ceptor,  including  the  encounter  time  and  the  relative 
motion  information.  Then,  the  maneuver  strategies  with 
different  conditions  must  be  much  different  regardless  of 
the  same  guidance  laws  and  maneuver  control  equations. 
Consequently,  we  have  to  design  the  respective  strategies 
when  the  encounter  time  is  determined  and  undeter¬ 
mined.  The  guidance  and  maneuver  logic  can  be  expressed 
in  Fig.  2. 

Fig.  2  shows  the  complete  guidance  and  maneuver  logic. 
To  begin  with,  as  the  basis  of  the  further  research  of  guidance 
law  and  maneuver  strategy,  it  needs  to  establish  the  3D 
nonlinear  coupled  relative  motion  equation  based  on  terminal 
constraints  and  current  motion  states.  Second,  the  nonlinear 
coupled  motions  equation  will  be  decoupled  into  linear 
ones  through  feedback  linearization,  and  the  3D  guidance 
law  which  can  satisfy  multiple  terminal  constraints  can  be 
also  obtained.  Furthermore,  it  is  easy  to  construct  maneuver 
control  equation  by  adding  maneuver  control  item  to  the 
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original  guidance  law.  With  this  equation,  it  will  establish  the 
optimal  performance  indices  to  evaluate  the  maneuver  per¬ 
formance  and  design  the  maneuver  strategies  using  optimal 
control  when  the  encounter  time  is  determined  and  unde¬ 
termined,  respectively.  In  addition,  considering  the  control 
capability  limitation,  especially  the  terminal  guidance  accu¬ 
racy,  the  performance  indices  under  different  situations  will 
be  reconstructed,  and  then  the  optimal  strategy  and  guaran¬ 
teed  cost  maneuver  strategy  will  be  obtained  as  well. 

This  section  will  take  full  advantage  of  optimal  control 
and  parameter  optimization  to  obtain  diving  optimal 
maneuver  strategy  with  the  determined  encounter  time. 

4.1.  Optimal  maneuver  strategy 


Reasonable  maneuver  strategy  is  necessary  for  the 
effective  penetration  [13],  The  penetration  performance 
is  closely  related  to  the  mutual  flight  velocities  and  over¬ 
loads  of  evader  and  pursuer  at  encounter  time.  As  the 
evader,  penetration  performance  increases  as  velocity  and 
overload  increase.  However,  it  is  always  unaccepted  to 
maneuver  blindly  because  of  too  much  energy  consump¬ 
tion  or  too  small  maneuver  overload  [13];  in  other  words, 
the  penetration  performance  must  be  declined  if 
maneuver  strategy  is  designed  unreasonable.  Conse¬ 
quently,  optimal  control  can  be  employed  to  resolve  this 
problem. 

To  begin  with,  take  longitudinal  maneuver  for  example; 
the  maneuver  control  equation  can  be  rewritten  as 
follows: 


J  x  =  —  kDx + UjnD 
\y=x 


(30) 


where  x  =  Z\  represents  the  longitudinal  LOS  angle  rate. 
We  can  introduce  optimal  performance  index  based  on 
energy  consumption  and  LOS  angle  rate  at  encounter  time: 

J  =  minJ-^2(t/)  +  2  £  “modtj  (31) 

where  tf  refers  to  the  encounter  time.  The  performance  index 
(31)  represents  the  minimum  value  of  the  sum  of  the 
maneuver  energy  consumption  in  the  time  interval  [0,  t/] 
and  the  maximum  value  of  LOS  angle  rate  at  encounter  time, 
which  means  we  should  design  a  control  input  umD  to  obtain 
maximum  velocity  and  overload.  With  the  performance  index 
(31),  the  optimal  control  law  of  Eq.  (30)  can  be  given: 
u*D  =  -K(t)x=  R  1  BrP(t)x  (32) 

in  which  B  =  [1],  F  =  [  - 1],  R  =  [1],  and  P(t)  can  be  calculated 
through  differential  Riccati  Eq.  [20]: 
j  P(t)  +  P(t)A+ArP(t)  -  P(t)BR  1  BrP(t)  =  0 
\P(.tf)  =  F  (  ) 


The  solution  of  Eq.  (33)  is 


exp(  -  2  kDtf)  +  (2  kD  -  l)exp(  -  2  kDt) 


(34) 


Different  encounter  times  tf  (5-10  s),  and  effective 
navigation  ratio  kD  =  4  are  set,  the  changing  curve  of 
optimal  gain  with  time  is  expressed  in  Fig.  3. 


t/(s) 


Fig.  3.  -K(t)  under  different  determined  encounter  times. 

Fig.  3  illustrates  that  the  maneuvering  flight  only  occurs 
near  the  encounter  time  to  save  energy  and  increase 
the  terminal  overload  as  much  as  possible.  Substituting 
Eq.  (34)  into  Eq.  (32),  the  optimal  maneuver  strategy  can 
be  computed: 

*  2kDexp(-2kDtf) 

mD  exp(  -  2kotf) +(2ko~1  )exp(  -  2kot) 


We  can  leam  from  the  optimal  guidance  gain  shown  in 
Fig.  3  and  Eq.  (35)  that  it  is  unnecessary  adopt  any  maneu¬ 
vering  flight  when  the  relative  distance  between  hypersonic 
vehicle  and  interceptor  is  large,  and  maneuver  only  occurs 
near  the  encounter  time  so  that  the  energy  can  be  saved  and 
terminal  overloads  can  be  increased  as  much  as  possible. 
However,  this  strategy  ignores  the  control  capability  limita¬ 
tion  and  the  influence  on  terminal  guidance  accuracy  caused 
by  maneuvering  flight,  which  means  that  the  maneuver 
performance,  control  capability  and  guidance  accuracy  can¬ 
not  be  considered  separately.  In  order  to  resolve  this  pro¬ 
blem,  we  have  to  reconstruct  the  performance  index  and 
strategy  based  on  the  above  analysis. 

The  performance  index  (31)  is  designed  to  get  mini¬ 
mum  energy  consumption  and  maximum  LOS  angle  rate  at 
encounter  time.  It  is  clear  that  minimum  energy  consump¬ 
tion  equals  to  maximum  terminal  velocity  magnitude,  and 
then  we  can  reconstruct  the  index  (31 )  as 

7=  sup  [(rfVEfx(t)+V®]  (36) 


In  Eq.  (36),  v0  is  the  initial  velocity:  the  velocity  is 
normalized  to  let  v(t)/vo  e  [0, 1],  kreViSe  is  the  revise  coeffi¬ 
cient  to  make  krevisex(t)  have  the  same  order  of  magnitude 
with  v(t)/v0.  The  performance  index  (36)  is  designed  with 
the  similar  meaning  with  index  (31 )  and  the  finite  time  Ha 0 
performance  index  [13],  but  it  will  be  more  effective 
because  practical  LOS  angle  rate  and  velocity  are  well 
considered.  It  is  clear  that  the  maneuver  performance 
will  be  optimal  when  weighted  average  of  the  velocity 
magnitude  and  LOS  angle  rate  at  encounter  time  is  the 
largest.  With  the  novel  performance  index  (36),  combining 
maneuver  performance,  control  capability  and  guidance 
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accuracy  together,  the  new  strategy  can  be  designed: 


where  I<Dmax  is  the  maximum  maneuver  amplitude,  g  can 
be  employed  to  modify  maneuver  start  and  finish  time. 
Being  similar  to  the  strategy  shown  in  Eq.  (35),  a  large 
amplitude  maneuver  near  the  encounter  time  is  most 
effective,  which  can  obtain  larger  overload  and  terminal 
velocity  magnitude.  Besides,  another  main  advantage  is 
that  the  maneuver  amplitude  will  reduce  to  zero  to  ensure 
the  terminal  guidance  accuracy,  and  the  rate  of  control 
gain  can  also  be  modified  by  modulating  the  maneuver 
parameters  KDmax  and  g.  Consequently,  compared  with  the 
strategy  shown  in  Eq.  (35),  the  maneuver  control  law  (37) 
is  more  usable  because  the  control  capability  and  terminal 
guidance  accuracy  are  fully  considered,  although  it  is  the 
suboptimal  maneuver  strategy. 

4.2.  Simulation  and  parameters  optimization 

In  the  simulation  of  guidance  performance  testing,  the 
1976  U.S.  Standard  Atmosphere  is  used  as  the  atmospheric 
model.  Lockheed-Martin's  CAV-H  vehicle  [1]  whose  weight 
is  907  kg,  reference  area  is  0.4839  m2  is  employed.  The 
equations  of  motion  are  integrated  with  a  fourth-order 
Runge-Kutta  method,  and  the  altitude  is  treated  as  the 
stopping  condition.  The  maneuver  control  parameters  are 
tf  =  60  s,  g  =  3,  and  krevise  =  180 Ik.  The  initial  and  terminal 
parameters  and  path  constraints  are  shown  in  Table  1. 

The  main  simulation  results  for  diving  guidance  and 
longitudinal  maneuver  are  shown  in  Figs.  4-7.  We  can 
learn  from  the  results  that  both  the  diving  high  precision 
and  maneuvering  flight  can  be  realized  with  the  aid  of 
proposed  strategy.  All  the  terminal  impact  angle  errors  are 
smaller  than  0.5°,  LOS  angle  rates  shrink  down  to  zero  at 
the  end  of  diving  flight,  and  additional  overloads  about  3- 
3.5  g  were  generated.  In  addition,  all  the  movement  states 
under  different  maximum  maneuver  amplitudes  keep  the 
same  when  the  fight  time  is  smaller  than  55  s  because  the 
amplitudes  are  zero  according  to  the  maneuver  strategy. 
When  fight  time  is  larger  than  55  s,  maneuvering  flight 
starts,  and  the  angles  of  attack  shown  in  Fig.  4  increase 


Table  1 

Simulation  parameters. 


Initial  velocity  magnitude  (v0) 

2000  m/s 

Initial  velocity  slope  angle  (S0) 

0° 

Initial  velocity  azimuth  (n0) 

45° 

Initial  longitude  (20) 

0° 

Initial  latitude  ( <pa ) 

0° 

Initial  altitude  (h0) 

40  km 

Terminal  longitude  [2f) 

1° 

Terminal  latitude  (q>f) 

1° 

Terminal  velocity  slope  angle  (0f) 

-80” 

Maximum  angle  of  attack  (amax) 

[  —  20°, 20°] 

Maximum  angle  of  attack  rate  (|«lmjx) 

5"  s-1 

Maximum  overload  (JVmax) 

[0, 20g] 

Maximum  bank  angle  rate  (lolmax) 

20' s-1 

t/(s) 


Fig.  6.  Longitudinal  overload  versus  time  curves. 


rapidly  to  generate  enough  aerodynamic  force  to  achieve 
large  amplitude  maneuver.  However,  as  shown  in  Fig.  5, 
due  to  the  maximum  angle  of  attack  and  the  rate  of  change 
limitation,  the  practical  LOS  angle  rates  and  overloads 
cannot  increase  as  required  amplitude  increases.  Con¬ 
sequently,  it  can  be  seen  from  Figs.  6  and  7  that  the 
difference  between  the  trajectories  with  different 
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maximum  amplitudes  is  not  so  distinct,  especially  when 
KDmax  >  0.2.  Moreover,  according  to  the  maneuver  strategy, 
the  amplitude  decreases  as  the  flight  time  increases  when 
t>  60  s  so  that  the  guidance  accuracy  can  be  ensured.  In 
general,  the  proposed  strategy  is  applicable  to  resolve  the 
diving  guidance  problem. 

It  is  known  that  maneuver  performance  can  be  influ¬ 
enced  by  maneuver  amplitude  and  start  time.  Therefore 
simulation  can  only  evaluate  that  the  new  proposed 
strategy  is  available  for  the  guidance  and  maneuver  tasks, 
but  the  performances  under  different  parameters  f  and 
KDmax  remain  to  be  clarified.  In  order  to  obtain  the  optimal 
parameters,  we  set  several  different  parameters  KDmax 
e  [0.1, 2.0],  f  e  [1,20]  under  the  same  initial  conditions  for 
each  case,  and  the  simulation  results  are  shown  in  Figs.  8 
and  9.  The  results  demonstrate  that  the  performance  index 
is  closely  related  to  f  and  KDmax',  both  of  them  cannot  be 
selected  too  small  or  too  large  because  of  too  much  energy 
consumption  or  too  small  LOS  angle  rate.  Combining  the 
results  in  Figs.  8  and  9,  it  can  be  concluded  that  the 
maneuver  performance  will  be  optimal  when  f  is  selected 
near  six.  In  addition,  or  each  given  £,  the  index  increases 
rapidly  as  the  parameter  KDmax  increases  when  it  is  smaller 
than  0.3.  However,  due  to  the  maximum  angle  of  attack 
and  the  rate  of  change,  the  maneuver  performance  cannot 
be  further  improved  when  Komax  is  larger  than  0.5.  As  a 
result,  we  can  obtain  elementary  conclusion  that  the 
maneuver  performance  will  reach  to  the  best  when  £  =  6, 
and  I<Dmax  is  larger  than  0.3. 

The  above  simulation  aims  to  achieve  maneuvering 
flight,  ignoring  the  influence  on  terminal  guidance  accu¬ 
racy.  The  indices  and  terminal  guidance  accuracies  under 
different  maneuver  parameters  f  and  I(Dmax,  including  miss 
distance  and  impact  angle,  are  shown  in  Table  2.  Firstly, 
compared  with  the  un-maneuverable  trajectory  with  the 
index  1.152,  the  maneuver  performance  can  be  much 
improved  by  maneuvering  flight.  Secondly,  setting  a  rela¬ 
tive  small  values  of  KDmax  and  f  can  achieve  a  highly 
precise  guidance,  otherwise  both  the  maneuver  perfor¬ 
mance  and  guidance  accuracy  will  be  degraded,  especially 
for  the  terminal  accuracy  Moreover,  as  a  complementarily 
for  the  above  simulation,  f=6  and  KDmax  =  0.4  are  the 


KDma 


Fig.  9.  Performance  indices  when  £  =  5, 6, 7. 


Guidance  accuracy  and  performance  index  under  different  f  and  Komax- 
KDmax  f  J  Miss  distance  (m)  Impact  angle  (deg) 


0 

0.2 

0.2 

0.2 

0.2 

0.4 

0.4 

0.4 

0.4 

0.6 

0.6 

0.6 

0.6 


1.152  1.652 

1.174  1.681 

1.187  1.124 

1.192  1.702 

1.185  2.221 

1.176  1.362 

1.189  1.743 

1.196  1.876 

1.194  56.21 

1.176  1.142 

1.190  1.798 

1.196  143.5 

1.194  416.7 


-80.097 

-79.873 

-79.933 

-80.203 


-79.731 

-80.809 

-82.330 

-79.772 

-80.126 

-86.413 

-87.979 


expected  optimal  maneuver  parameters  owing  to  the 
dramatic  increase  of  miss  distance  and  impact  angle  error 
increase  dramatically  with  a  I(Dmax  larger  than  0.6. 
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5.  Optimal  maneuver  strategy  when  encounter  time  is 
undetermined 

However,  as  for  the  hypersonic  vehicle  which  is 
regarded  as  an  evader,  the  design  of  maneuver  strategy 
will  be  much  more  difficult  if  interceptor's  movement 
information,  such  as  the  encounter  time,  cannot  be 
determined  perfectly.  Guaranteed  cost  maneuver  strategy 
which  needs  maneuvering  flight  to  be  carried  out  all  over 
the  dive  phase  can  be  adopted  to  resolve  this  problem.  Be 
similar  to  the  above  optimal  strategy,  optimal  control  will 
be  employed  to  design  the  maneuver  strategy  here,  term¬ 
inal  guidance  accuracy  and  control  capability  will  be  fully 
considered,  and  the  corresponding  parameters  will  be  also 
optimized. 

5.1.  Guaranteed  cost  maneuver  strategy 


Tg/(s) 


Fig.  10.  -K (t)  under  different  terminal  times. 


Taking  longitudinal  maneuver  for  example,  in  order  to 
obtain  the  guaranteed  cost  maneuver  strategy  based  on 
optimal  control,  we  should  introduce  optimal  performance 
index  at  first: 

J  =  min  ^  x(Tg)Fx(Tg) + \  ‘  -  x2  dtj  (38) 

where  Tg  is  the  flight  time-to-go.  The  performance  index 
(38)  is  constructed  to  obtain  minimum  energy  consump¬ 
tion  and  maximum  LOS  angle  rate  during  the  whole  dive 
phase,  and  let  terminal  LOS  angle  rate  equal  to  the 
required  one.  With  the  state  space  Eq.  (30)  and  perfor¬ 
mance  index  (38),  the  key  step  of  obtaining  optimal 
control  law  is  the  solution  of  differential  Riccati  Eq.  [20]: 

P(t)+P(t)A+ATP(t)-P(t)BR-'BTP(t)+<l(t)  =  0 

P(tf)  =  F  (39) 

where  A  =  [  -  kD],  B = ft],  R  =  [1],  Q  =  [- 1],  F»  [0],  the  solu¬ 
tion  of  Eq.  (39)  is 

(Tg  - 1)  1  +  atanh  j  j 

(40) 

substituting  P(t)  into  optimal  control  law,  the  maneuver 
strategy  can  be  obtained. 


P(t)  =  -  kD  +  \/k$)- 1  ■  tanh 


u*  D  =  -  I<(t)x  =  ^fcD  -  \JkA0-l  ■  tanh 

(Tg  - 1) \Jk2D-l  +  atanh  ^  — j  j  j  x  (41) 

Different  guidance  finish  time  Tg(50-100  s),  and  effective 
navigation  ratio  kD  =  4  are  set;  the  changing  curve  with 
time  is  expressed  in  Fig.  10. 

As  shown  in  Fig.  10,  we  can  leam  that  the  maneuvering 
flight  should  be  carried  out  during  the  whole  flight  process  to 
ensure  the  maneuver  performance,  and  maneuver  control 
gain  will  reduce  to  zero  rapidly  to  decrease  the  influence  on 
terminal  guidance  accuracy  as  much  as  possible.  However, 
this  strategy  can  only  improve  maneuver  performance  by 


increasing  the  terminal  velocity  instead  of  LOS  angle  rate 
because  of  the  constant  maneuver  gain,  ignoring  the  control 
capability  limitation  and  the  path  constraints  resulting  from 
the  too  large  rate  of  the  control  gain  in  the  end  of  dive  phase. 
Consequently,  we  have  to  reconstruct  the  performance  index 
and  the  theoretical  optimal  maneuver  strategy. 

Maneuver  performance  is  closely  related  to  the  flight 
velocity  and  overload,  which  means  that  with  the  larger 
the  velocity  and  overload,  the  better  the  maneuver  per¬ 
formance.  Therefore,  for  similar  performance  index  (38), 
we  can  reconstruct  the  performance  index  which  consists 
of  the  terminal  velocity  magnitude  and  the  LOS  angle  rate 
integration  during  the  whole  dive  phase  to  make  it  more 
usable  in  practice.  Here,  the  LOS  angle  rate  integration  is 
employed  to  describe  the  overload  magnitude. 

Jr  =  sup  + k^  j  ifldtj  (42) 

in  which  v0  and  v(Tg)  are  the  initial  and  practical  terminal 
velocities,  respectively,  and  kv,kj  are  the  revise  coeffi¬ 
cients.  This  performance  is  designed  to  obtain  maximum 
terminal  velocity  and  overload  which  can  be  increased  by 
maneuvering  flight.  However,  it  is  known  that  maneuver¬ 
ing  flight  can  cause  serious  influence  to  terminal  guidance 
accuracy,  and  even  lose  the  target  if  the  amplitude 
is  designed  too  large.  Consequently,  based  on  the  novel 
constructed  index  (42),  it  is  better  to  construct  the 
integrated  performance  index  combining  terminal  gui¬ 
dance  accuracy  together. 

J2  =  sup  [k^^*hJo  *Apdt-k£rEr-k9|0(Tg)-6>/|J  (43) 

where  Er  represents  terminal  miss  distance,  0(Tg)  is  the 
practical  terminal  velocity  slope  angle,  both  kEr  and  k„ 
are  the  corresponding  revise  coefficients.  Compared  with 
the  index  which  only  aims  to  obtain  minimum  energy 
consumption  and  maximum  predict  error  described  in 
Ref.  [13],  this  index  is  more  usable  because  both  the 
maneuver  performance  and  guidance  accuracy  are  fully 
considered.  This  index  shows  that  the  larger  the  terminal 
velocity  and  XD  dt,  the  better  the  integrated  perfor¬ 
mance.  By  contrast,  the  smaller  the  terminal  impact  errors, 
the  better  the  integrated  performance.  As  a  result,  we  have 
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to  design  the  maneuver  strategy  to  improve  maneuver 
performance  and  ensure  guidance  accuracy. 

In  order  to  increase  the  reconstructed  performance 
index  (43)  as  much  as  possible,  it  is  necessary  to  recon¬ 
struct  the  maneuver  strategy  further.  Great  many  maneu¬ 
ver  models  can  be  well  employed  to  enhance  maneuver 
performance,  such  as  constant  acceleration  model,  Singer 
model  and  sinusoid  model  which  can  be  also  called  weave 
model,  in  which  the  last  one  is  the  most  effective  com¬ 
monly  [18],  Weaving  maneuver  has  already  been  adopted 
to  improve  maneuver  performance  for  anti-ship  missile 
successfully  [21],  but  it  cannot  be  employed  to  achieve 
maneuver  for  hypersonic  vehicle  in  dive  phase  directly 
because  of  the  serious  influence  on  the  terminal  guidance 
accuracy.  Therefore,  in  order  to  take  full  advantage  of  the 
excellent  performance  of  weaving  maneuver  and  decrease 
the  influence  on  the  guidance  accuracy  at  the  same  time, 
which  means  generate  flexible  trajectory  and  enhance 
maneuver  performance,  we  should  redesign  the  maneuver 
strategy  based  on  the  weaving  trajectory  proposed  in 
Refs.  [18,21], 


( 


h  >  hmidd 
h  <  hmidd 


(44) 


where  KDmax  and  FD  are  the  maneuver  amplitude  and 
frequency,  respectively,  h  is  the  practical  altitude  and  hmidd 
is  the  finish  altitude  of  maneuver.  The  reconstructed 
maneuver  strategy  can  be  explained  as  weaving  maneuver 
will  be  adopted  to  enhance  the  maneuver  performance 
when  h  >  hmidi,  and  the  maneuver  amplitude  keeps  con¬ 
stant  zero  if  h<hmidd  to  ensure  the  terminal  guidance 
accuracy.  In  this  way,  the  maneuver  performance  can  be 
improved  with  the  weave  maneuver  model,  and  the 
guidance  accuracy  can  also  be  ensured  because  of  the  zero 
maneuver  amplitude  when  h<hmidd.  Like  the  optimal 
strategy  with  determined  encounter  time,  this  strategy  is 
the  suboptimal  solution  of  the  guaranteed  cost  maneuver 
strategy,  but  it  will  be  more  useful  because  the  terminal 
guidance  accuracy  is  also  fully  considered.  Besides,  differ¬ 
ent  maneuver  trajectories  can  also  be  obtained  through 
modifying  the  control  parameters. 


5.2.  Simulation  and  parameters  optimization 

The  same  initial  and  terminal  simulation  parameters 
introduced  in  Section  4.2  are  taken,  and  CAV-H  model  is  also 
employed  in  the  simulation.  The  main  simulation  results  with 
the  maneuver  parameters  FD  =  4,  hmidd  =  5  km,  kv  =  1, 
fq  =  1 80 /n,kEr  =  0.01  and  k9  =  1  are  shown  in  Figs.  11-14.  It 
can  be  seen  from  the  results  that  the  novel  proposed  3D 
guidance  law  and  guaranteed  cost  maneuver  strategy  can 
steer  the  vehicle  to  attack  the  target  with  the  required  angle 
even  if  large  amplitude  maneuvering  flights  are  adopted  when 
hmidi  >  5  km-  The  terminal  miss  distances  and  impact  angle 
errors  are  smaller  than  2  m  and  1°,  respectively.  The  maneu¬ 
ver  trajectories  including  angle  of  attacks  are  much  different 
from  the  trajectories  without  maneuver  and  they  will  be  more 
flexible  when  the  amplitude  gets  larger.  During  the  whole 
dive  phase,  there  are  great  many  collapses  for  the  LOS  angle 
rates  shown  in  Fig.  12,  overloads  expressed  in  Fig.  13  and 


velocity  slope  angles  depicted  in  Fig.  14,  but  the  terminal 
guidance  accuracy  can  be  ensured  because  of  the  zero 
maneuver  amplitude  at  the  end  of  flight  However,  due  to 
the  maximum  angle  of  attack  and  the  limited  rate  of  change, 
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Fig.  15.  Maneuver  performance  indices. 


the  practical  maneuver  ranges  cannot  increase  as  the  man¬ 
euver  amplitudes  increase.  These  simulation  results  indicate 
the  practicability  of  the  3D  guidance  law  and  the  proposed 
guaranteed  cost  maneuver  strategy. 

The  above  simulation  can  only  confirm  the  guaranteed 
cost  maneuver  strategy  can  accomplish  high  precision  gui¬ 
dance  as  well  as  maneuvering  flight  at  the  same  time,  but 
the  maneuver  performance  under  different  parameters  still 
needs  to  be  validated.  In  order  to  test  the  maneuver 
performance  under  different  maneuver  frequencies  and 
amplitudes  and  obtain  the  optimal  guaranteed  cost  maneu¬ 
ver  strategy,  we  set  several  different  maneuver  parameters 
under  the  same  initial  condition,  including  KDmax= 
0-1  and  Fd= 0-5.  The  maneuver  performance  indices  and 
integrated  performance  indices  are  shown  in  Figs.  15  and  16, 
respectively. 

Fig.  15  shows  that  the  maneuver  performance  which  is 
closely  related  to  the  amplitude  and  frequency  which  can 
be  enhanced  by  maneuvering  flight  effectively.  For  each 
given  frequency  (KD  max  t4  0),  the  amplitude  of  maneuver 
performance  increases  initially  and  decreases  afterwards 
to  a  nearly  same  value  due  to  the  large  energy  consump¬ 
tion.  The  peak  values  of  maneuver  performance  are  much 
different  from  each  other  for  different  given  frequencies. 
Additionally,  it  indicates  that  the  smaller  the  maneuver 
frequency  is,  the  better  the  maneuver  performance  can  be 
achieved  when  the  amplitude  is  less  than  0.6.  The  main 
reason  is  that  the  maneuver  range  and  the  LOS  angle  rate 
in  each  maneuver  period  increase  as  the  frequency 
decreases  when  the  amplitude  is  small.  Flowever,  the 
performance  with  different  frequencies  tends  to  decrease 
when  the  amplitude  is  above  0.6  because  of  the  large 
energy  consumption  and  control  capability  limitation. 
Consequently,  reasonable  maneuver  frequency  and  ampli¬ 
tude  are  necessary  to  improve  the  maneuver  performance. 

Fig.  16  and  Table  3  illustrate  that  the  integrated  perfor¬ 
mance  indices  are  much  different  when  the  terminal 
guidance  accuracy  is  fully  considered,  and  maneuvering 
flight  can  exert  a  serious  influence  on  the  terminal 
guidance  accuracy,  especially  for  large  amplitude  maneu¬ 
ver.  Owing  to  the  improvement  of  maneuver  performance 
and  relative  high  guidance  accuracy,  the  integrated  per¬ 
formance  index  under  each  given  maneuver  frequency 


Fig.  16.  Integrated  performance  indices. 


Guidance  accuracy  and  performance  index  under  different  KDwax  and  FD. 


KDmax  Fd  J,  J2  Miss  Impact 

distance  (m)  angle  (deg) 


0  -  2.529 

0.2  2  1.858 

0.2  4  1.985 

0.2  6  1.992 

0.4  2  11.34 

0.4  4  2.213 

0.4  6  4.201 

0.6  2  8.131 

0.6  4  17.15 

0.6  6  2.916 

0.8  2  7.690 

0.8  4  11.700 

0.8  6  12.297 


2.429  1.648 

1.561  1.062 

-5.494  283.2 

1.996  1.428 

4.159  1.282 

-37.203  403.7 

-30.15  98.66 

2.849  1.240 

-45.824  797.6 

-3.700  31.74 

-1.063  38.07 


-80.096 

-80.295 

-80.391 

-80.009 

-88.273 

-80.216 

-79.959 

-83.127 

79.904 

-79.936 

-81.919 

-84.593 

-84.485 


(Kd  max  #  0)  increases  initially  with  increasing  amplitude. 
But  the  index  turns  to  decrease  if  the  amplitude  continues 
increasing  because  of  the  decline  of  guidance  accuracy 
caused  by  large  amplitude  maneuver.  Moreover,  for  the 
given  amplitude,  the  larger  the  maneuver  frequency  is,  the 
better  the  integrated  performance  index  can  be  achieved, 
especially  for  relatively  small  amplitude.  The  main  reason 
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is  that  the  maneuver  range  and  the  errors  in  each  period 
caused  by  maneuvering  flight,  including  position  and 
direction  errors,  will  be  enlarged  when  the  frequency  is 
small.  It  is  these  errors  which  cannot  be  eliminated  by  the 
limited  control  capability  that  leads  to  the  decline  of  the 
terminal  guidance  precision  and  performance  index.  In 
addition,  the  guidance  accuracy  and  integrated  index  will 
continue  declining  as  the  amplitude  increases. 

The  simulation  results  for  the  two  different  performance 
indices  are  compared;  we  can  leam  that  contradiction  exists 
between  high  precision  guidance  and  maneuvering  flight. 
For  example,  the  maneuver  performance  is  the  optimal 
when  I<Dmax  =  0.6,  FD  =  4,  but  integrated  performance  is 
decreased  so  much  because  of  the  too  large  impact  errors. 
Therefore,  both  the  maneuver  frequency  and  amplitude 
cannot  be  selected  too  large  or  too  small  no  matter  whether 
terminal  guidance  accuracy  is  considered. 

6.  Conclusions 

This  study  presents  a  novel  diving  maneuver  strategy 
considering  terminal  guidance  accuracy  for  hypersonic 
vehicles.  The  main  characteristic  of  the  proposed  metho¬ 
dology  lies  in  the  fact  that  it  takes  full  advantage  of  three- 
dimensional  diving  movement,  optimal  control  and  para¬ 
meter  optimization  to  realize  high  precision  guidance  and 
obtain  optimal  maneuver  strategy  at  the  same  time.  In 
addition,  considering  the  control  capability  limitation,  the 
performance  indices  and  maneuver  strategies  are  recon¬ 
structed  when  the  encounter  time  is  determined  and 
undetermined.  The  major  advantages  of  the  proposed 
algorithm  can  be  listed  as  follows: 

1 )  The  three-dimensional  nonlinear  relative  coupled  motion 
equations  are  constructed  without  any  assumptions  and 
approximations,  which  means  the  coupling  characteristic 
is  fully  materialized  in  the  equations. 

2)  Both  the  feedback  linearization  and  feedback  control  are 
adapted  so  that  the  solution  of  the  guidance  problem  is 
greatly  simplified,  and  the  maneuver  control  equations 
are  obtained  based  on  the  diving  guidance  law. 

3)  It  makes  best  use  of  optimal  control  to  obtain  the 
maneuver  strategies,  and  reconstructs  the  performance 
indices  and  strategies  considering  terminal  guidance 
accuracy  and  control  capability.  All  the  maneuver  control 
parameters  are  optimized  by  simulation  experiment. 

4)  It  can  be  concluded  that  the  maneuvering  flight  near 
the  encounter  time  is  the  most  effective  strategy.  The 
guaranteed  cost  maneuver  strategy  is  necessary  when 
the  encounter  time  is  undetermined,  and  both  the 
maneuver  frequency  and  amplitude  cannot  be  selected 
too  small  or  large. 

However,  there  are  also  some  disadvantages,  in  which 
the  maneuver  control  parameters  optimization  is  the  main 


research  topic  for  further  study.  In  general,  the  proposed 
algorithm  can  provide  some  reference  for  the  diving 
guidance  of  an  RLV,  hypersonic  vehicle,  etc. 
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